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Abstract: Effects of mowing on the composition and diversity of grasslands varied with climate change 
(e.g., precipitation and temperature). However, the interactive effects of long-term mowing and climate 
change on the diversity and stability of leguminous and non-leguminous species in the semi-arid 
grasslands are largely unknown. Here, we used i situ monitoring data from 1982 to 2011 to examine the 
effects of continuous mowing and climate change on the plant biomass and diversity of leguminous and 
non-leguminous species, and soil total nitrogen in the typical semi-arid grasslands of northern China. 
Results showed that the biomass and diversity of leguminous species significantly decreased with the 
increasing in the biomass and diversity of non-leguminous species during the 30-a period. Variations in 
biomass were mainly affected by the long-term mowing, while variations in diversity were mainly explained 
by the climate change. Moreover, the normalized change rates of diversity in leguminous species were 
significantly higher than those in non-leguminous species. Mowing and temperature together contributed 
to the diversity changes of leguminous species, with mowing accounting for 50.0% and temperature 
28.0%. Temporal stability of leguminous species was substantially lower than that of non-leguminous 
species. Consequently, soil total nitrogen decreased in the 2000s compared with the 1980s. These findings 
demonstrated that leguminous species were more sensitive to the long-term mowing and climate change 
than non-leguminous species in the semi-arid grasslands. Thus, reseeding appropriate leguminous plants 
when mowing in the semi-arid grasslands may be a better strategy to improve nitrogen levels of grassland 
ecosystems and maintain ecosystem biodiversity. 
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1 Introduction 


Rapid loss of global biodiversity in grasslands, the largest terrestrial ecosystem worldwide, has 
greatly affected community composition and ecosystem functioning at a broad scale (Tilman et al., 
2006; Cardinale et al., 2012; Hautier et al., 2015). Previous experimental studies conducted in 
extensively natural grassland systems have demonstrated that ongoing decrease in biomass 
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production is associated with a decline in plant species diversity (Tilman, 1999; Roscher et al., 
2005; Cardinale et al., 2007). Explanations for this trend focus on hypotheses that less diverse 
communities occupy smaller niche spaces. Therefore, the community cannot fully utilize 
distributed resources, thus reducing plant—plant interactions (He et al., 2016). Consisting with 
these hypotheses, a loss of biodiversity also indicates that important plant species or plant 
functional groups have been lost (Cardinale et al., 2012; Chen et al., 2016). 

A large body of research has examined climatic and anthropogenic factors as causal 
mechanisms underlying diversity loss and productivity declines in the grasslands (Moritz and 
Agudo, 2013; Hautier et al., 2015). For example, simulated precipitation experiments have shown 
that precipitation has an important effect on grass and forb species diversity and ecosystem 
production (Baez et al., 2013; Gherardi and Sala, 2015). Nevertheless, a study concluded that 
temperature was a better indicator than precipitation as a predictor of plant community assembly 
across dryland regions (Butterfield and Munson, 2016). Lin et al. (2010) found that warming 
significantly increased biomass by 12.3% across all the terrestrial plants. This terrestrial plant 
biomass response to warming would potentially lead to changes in community composition and 
structure, biodiversity and ecosystem functioning in a warmer environment. In addition to 
climatic factors, anthropogenic activities may play more important roles in altering grassland 
diversity and productivity relationships (Zhou et al., 2006; Miko and Storch, 2015). As one of the 
major management techniques for natural grasslands, mowing has markedly altered the 
composition and diversity of vegetation communities (Yang et al., 2012; Talle et al., 2016; Liu et 
al., 2018). Elucidating the effects of mowing, climate change and their possible interactions on 
productivity and diversity will help conserve biodiversity and stability of grassland ecosystems 
under global change scenarios. 

Vegetation communities consist of various functional groups, including legumes and 
non-N>-fixing grasses and forbs (Cowles et al., 2016; Korell et al., 2016). Leguminous species, a 
crucial plant functional group in grassland ecosystems, are critical for maintaining grassland 
sustainability and performing ecosystem functions (e.g., biomass production and nutrient cycling) 
(Wu et al., 2016; Yang et al., 2016; Gao et al., 2017). Meanwhile, nitrogen, which is symbiotically 
fixed by legumes, is one of the key drivers of the positive relationship between productivity and 
biodiversity (Li et al., 2016). Yang et al. (2011) found that climate change, e.g., increased 
precipitation, could differentially affect the abundance of legumes and non-legumes in a 6-a 
studied period. However, to date, there has been a lack of studies examining the response of 
leguminous plant functional groups to human disturbances (e.g., long-term mowing and climate 
change), posing serious challenges for predicting vegetation community dynamics under 
projected global climate change scenarios. 

The temperate steppe of northern China, a typical biome consisting of semi-arid regions, has an 
important and long-term provision capacity for people living, and it has experienced a great deal 
of anthropogenic disturbances and climate change (Lavorel et al., 1997; Peng et al., 2013). 
However, the effects of these changes on plant species sensitivity between long-term and 
short-term studies are different. For example, long-term mowing increased plant diversity (Yang 
et al., 2012), but short-term mowing did not affect vegetation community structure in a semi-arid 
steppe (Liu et al., 2018). Here, our study was conducted to evaluate the effects of long-term 
mowing and climate change on functional group sensitivity of leguminous and non-leguminous 
species using detailed monitoring data from 1982 to 2011. Specifically, the following two 
questions will be focused: (i) why do leguminous and non-leguminous species respond to climate 
change and mowing in a different way? and (ii) which functional group respond stronger to the 
environmental change? 


2 Materials and methods 


2.1 Study area 
This study was conducted at the Inner Mongolia Grassland Ecosystem Research Station 
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(IMGERS; 43°38'N, 116°42’E; 1200 m a.s.l.) of the Chinese Academy of Sciences, located in the 
Xilin Gol League of Inner Mongolia Autonomous Region, China. This region is characterized by 
a temperate semi-arid continental climate, with an annual mean temperature of approximately 
0.5°C and a mean annual (1982-2011) precipitation of approximately 350 mm, 60%-80% of 
which falls during the growing season (June-August; Baoyin et al., 2014) and coincides with the 
highest annual temperatures (Chen et al., 2016). The site has a sandy-loam dark chestnut soil 
(Calcic Chernozem according to the US soil taxonomy classification system), including a humus 
layer of 20-30 cm and a calcic layer at approximately 50 cm in depth (Bai et al., 2004). Before the 
experiment began, the grasslands had experienced mowing for haymaking and grazing by sheep 
and cattle at a low stocking rate. The dominated plant was Leymus chinensis and comprised 55% 
(+15%) of the total hay production over the years, with Stipa grandis and legumes comprising 
25% (412%) and 5% (+4%), respectively (Baoyin et al., 2014). 


2.2 Mowing experiment 


Mowing experiment conducted from 1982 to 2011 was previously described in detail by Baoyin 
et al. (2014). The mowing experiment consisted of four treatments, i.e., mowing once every other 
year; mowing once per year for two successive years, with no mowing for the third year; mowing 
once every year and mowing twice every year that have been continually maintained since 1982. 
No mowing was considered as the control treatment. The long-term monitoring data of two 
treatments were used in this study, i.e., mowing once every year and no mowing treatment. Each 
treatment had six replicates that were randomly assigned to six plots. Each plot was 2 mx2 m, and 
there were 12 plots in total. Plots were also separated by 1-m walkways. 


2.3 Variables for community structure 


All plants were cut and kept to 6 cm height aboveground over each plot at each mowing, but only 
the plants at the center of each plot were sampled to determine the aboveground plant biomass 
and species diversity. Plants were sorted into leguminous and non-leguminous species, oven-dried 
at 65°C for 48 h and weighed. Plants were approximately mowed on 16 August from 1982 to 
2011 (the specific mowing date was determined by the meteorological conditions of the current 
year) (Baoyin et al., 2014). 


2.4 Soil sampling and analysis 


In the 1980s and 2000s, three topsoil cores (3 cm in diameter and 30 cm in depth) were randomly 
collected from each plot and then mixed to form one uniform soil sample per plot. After the soil 
was gently mixed and roots were removed, the moist soil was passed through a 2-mm mesh sieve, 
air-dried and used to determine the total nitrogen (TN) with a Kjeltec analyzer (Kjeltec 2300 
Analyzer Unit, Foss, Sweden). All results are expressed on a dry matter basis. 


2.5 Meteorological data collection 


From 1982 to 2011, meteorological factors including monthly mean temperature and mean 
monthly precipitation at the experimental site were recorded by IMGERS. The dynamic 
characteristics of temperature and precipitation during the growing season from 1982 to 2011 are 
shown in Figure S1. 


2.6 Statistical analysis 


Aboveground biomass (AGB), relative aboveground biomass (RAGB), Shannon-Wiener diversity 
index (H), and relative Shannon-Wiener diversity index (RH) of leguminous and non-leguminous 
species were calculated. We calculated the H as follows: H= —)‘p;xIn(pi), where p; is the relative 
abundance of species i. Specifically, RAGB and RH are the proportions of AGB and H of 
leguminous or non-leguminous species in the total biomass or diversity of the whole community, 
respectively. The diversity effect is the ratio of AGB to H in leguminous and non-leguminous 
species. Change rate (CR) of the diversity effect was calculated by the slope between the diversity 
effect and year. 

A one-way ANOVA (analysis of variance) with Duncan's multiple-range tests was performed to 
test the effect of mowing on diversity from 1982 to 2011. We also used a linear regression to 
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determine the dynamics of AGB, RAGB, H and RH of leguminous and non-leguminous species. 

In addition, the AGB, RAGB, H and RH values of leguminous and non-leguminous species 
were further normalized, and normalized change rate (NCR) was used to assess the sensitivity of 
leguminous and non-leguminous species to mowing. The normalized values were calculated using 
the following equation (Karhu et al., 2014): 


DN D igin ~ Dsn (1) 

7 D nean ? 
where DN, Dorigin, and Dmean are the normalized, original and mean values of AGB, RAGB, H and 
RH of leguminous and non-leguminous species, respectively. Further, we calculated NCR values 
of AGB, RAGB, H and RH of leguminous and non-leguminous species from 1982 to 2011 using 


the least squares method given by the following equation (Karhu et al., 2014): 
NCR = Sow, (2) 


year 
where Apn and Ayear are the changed normalized values and changed year values of AGB, RAGB, 
H and RH of leguminous species and non-leguminous species, respectively. 
We analyzed the degrees of temporal stability (TS) of leguminous and non-leguminous 
species in the sampled communities using a stability index that was represented by the reciprocal 
of the coefficient of variation (CV; Yang et al., 2017). 


Tso (3) 


where SD and Dmean are the standard deviation and mean values of AGB in the main species as 
well as the functional groups of leguminous and non-leguminous species, respectively. 

Structural equation modeling (SEM) was performed to analyze different hypothetical pathways 
that may explain the linkages between mowing and climatic factors (including temperature and 
precipitation) to leguminous and non-leguminous species. In the SEM approach, we developed 
preliminary models based on hypothesized relationships between AGB and diversity of 
leguminous and non-leguminous species with mowing and climatic factors (Li et al., 2015). Then, 
according to the significant regression weights, it was also necessary to change the relationship 
among variables. The utility of AGB and diversity of leguminous and non-leguminous species 
within the SEM was compared through a number of measures. The maximum likelihood y? 
goodness-of-fit tests and Akaike information criterion (AIC) were used to examine the 
goodness-of-fit of the preliminary modeling. 

We conducted all standard statistical analysis using SPSS software version 19.0 (IBM Corp., 
Armonk, NY, USA) at P<0.05 and P<0.01 levels. The SEM analysis was performed by the 
AMOS 18.0 software package (SPSS Inc., Chicago, IL, USA). All graphs were prepared using 
SigmaPlot version 12.0 (Systat Software Inc., San Jose, CA, USA). 


3 Results 


3.1 AGB and species diversity 


The AGB and RAGB of leguminous species gradually declined from the 1980s to the 2010s 
(Fig. 1a and b). Mowing significantly decreased AGB by 54.43% (P<0.01) and RAGB by 25.26% 
(P<0.05) of leguminous species compared with control (Fig. 1a and b). In addition, mowing 
decreased AGB but significantly increased RAGB of non-leguminous species (Fig. S2). 
According to normalized analyses, the absolute biomass and relative biomass of NCR of 
leguminous species were significantly higher than those of non-leguminous species (Fig. 1c and 


d). 

The significant decreases in H and RH of leguminous species were also observed from 1982 to 
2011 (Fig. 2a and b), while an increase in RH of non-leguminous species was found (Fig. S2). 
The absolute biomass and relative biomass of NCR of leguminous species were significantly 
higher than those of non-leguminous species under mowing, but there were no significant 
differences between them under control (Fig. 2c and d). 
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Fig. 1 Dynamics of biomass (a) and relative biomass (b) of leguminous species from 1982 to 2011, and 
normalized change rate (NCR) values of absolute biomass (c) and relative biomass (d) of leguminous and 
non-leguminous species under mowing and control. Relative biomass of leguminous species is the proportion of 
AGB of leguminous species to the total biomass of community. Bars are standard errors. Different lowercase 
letters represent significant differences between leguminous and non-leguminous species at P<0.05 level. 
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Fig. 2 Dynamics of species diversity (a) and relative species diversity (b) of leguminous species from 1982 to 
2011, and normalized change rates (NCR) of absolute diversity (c) and relative diversity (d) of leguminous and 
non-leguminous species under mowing and control. Relative species diversity of leguminous species is the 
proportion of species diversity of leguminous species to the total species diversity of community. Bars are 
Different lowercase letters indicate significant differences between leguminous and 
non-leguminous species at P<0.05 level. 
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SEM analysis was further used to determine the factors driving the sensitivity of leguminous species. 
Variations in AGB could be explained by mowing, with 39.0% of the total variation for leguminous 
species and 52.0% for non-leguminous species (Fig. 3a and b). However, variations in species 
diversity could be explained by climatic factors, with 56.0% of the total variation for leguminous 
species and 67.0% for non-leguminous species (Fig. 3 a and b). For the relative biomass, about 70.0% 
of variation could be indirectly explained for leguminous species (Fig. 3c), with climatic factors 
accounting for 63.4% and mowing 7.5% (Fig. 3d). 
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Fig. 3 Pathways for the analysis of changes in leguminous (a) and non-leguminous species (b) biomass, relative 
biomass decline (c) and net effect value on relative biomass decline of leguminous species (d) under climate 
change and mowing. Values with arrows represent standardized path coefficients. Solid and dashed arrows 
represent significant and non-significant paths in a fitted structural equation model (a, y7=0.96, P=0.81; b, 77=2.304, 
P=0.52; c, v7=20.28, P=0.11). ***, P<0.001 level; **, P<0.01 level; *, P<0.05 level; ns, P>0.05 level; PR, 
precipitation; TE, temperature; CF, climatic factors; M, mowing; LEp, species diversity of leguminous species; 
LEm, biomass of leguminous species; NLp, species diversity of non-leguminous species; NLm, biomass of 
non-leguminous species; LRw, relative biomass decline in leguminous species. 


3.2 Relationship between AGB and H 


There was a significantly positive correlation between AGB and H of leguminous and non-leguminous 
species under mowing (Fig. 4). The H of leguminous species significantly decreased over time, but 
there was no change for non-leguminous species (Fig. 5a and b). Mowing significantly decreased H by 
68.28% compared with control (Fig. 5a). The CR of biomass/diversity ratio under mowing for 
leguminous species was much higher than that for non-leguminous species under control (Fig. 5c). The 
NCR of biomass/diversity ratio for leguminous species was significantly higher than that for 
non-leguminous species under both mowing and control (Fig. 5d). Furthermore, the divergent effects 
of H of leguminous and non-leguminous species could be explained by the differences in their 
responses to mowing and climatic factors (Fig. 6; Table 1). The H of leguminous species was mainly 
explained by both mowing and temperature, with mowing accounting for 50.0% and temperature 
28.0% (Fig. 6a). However, the H of non-leguminous species was mainly explained by mowing, with 
54.0% of the total variation (Fig. 6b). 


202011.00135v1 


chinaXiv 


ChinaXivA ERAT! 


JOURNAL OF ARID LAND 


) 


șa 300 


R°=0.18, P<0.05 
R°=0.19, P<0.05 


@ Control 
@ Mowing 


R*=0.39, P<0.01 


non-leguminous species (g/m 


Aboveground biomass of 
leguminous species (g/m 


Aboveground biomass of 


0 
4 6 8 10 12 14 16 18 20 22 
Species diversity of leguminous species Species diversity of non-leguminous species 


Fig. 4 Relationships between aboveground biomass and species diversity of leguminous (a) and non-leguminous 
(b) species 
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Fig. 5 Biomass/diversity ratios of leguminous (a) and non-leguminous species (b) from 1982 to 2011. And 
change rate (CRbiomass/diversity, C) and normalized change rate (NCRobiomass/diversity, d) of biomass/diversity ratios of 
leguminous and non-leguminous species under mowing and control. Bars are standard errors. Different lowercase 
letters represent significant differences between leguminous and non-leguminous species at P<0.05 level. 


Fig. 6 Structural equation modeling analysis of the pathways representing changes in the diversity effects of 
leguminous (a) and non-leguminous species (b). Values with arrows represent standardized path coefficients. 
Solid and dashed arrows respectively represent significant and non-significant paths in a fitted structural equation 
model (a, y*=0.01, P=0.99; b, 7°=0.01, P=0.99). ***, P<0.001 level; *, P<0.05 level; ns, P>0.05 level. PR, 
precipitation; TE, temperature; M, mowing; DEL, diversity effects of leguminous species; DEN, diversity effects 
of non-leguminous species. 
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Table 1 Regression coefficients of biomass and species diversity of leguminous species with climatic factors 
under control and mowing treatments from 1982 to 2011 


Precipitation Temperature 
Parameter Treatment 
Linear Nonlinear Linear Nonlinear 

Control 0.05 0.34 —0.50" —0.50* 

Biomass m a 
Mowing 0.14 0.30 —0.59 —0.60 
Control —0.03 0.42 —0.32 -0.33 

Relative biomass a n . 
Mowing 0.02 0.54 0.57 —0.50 

Control 0.25 0.31 —0.56™ -0.58" 

Species diversity , . 
Mowing 0.09 0.44 —0.46 —0.64 

Control 0.22 0.19 -0.39" -0.62™" 

Relative species diversity : E 
Mowing 0.03 0.32 —0.38 —0.60 


Note: **, P<0.01 level; *, P<0.05 level. 


3.3 Temporal stability 


Mowing slightly decreased the temporal stability of AGB of leguminous and non-leguminous species 
(Fig. 7a). But the temporal stability of AGB in leguminous species was substantially lower than that of 
non-leguminous species (Fig. 7a). There was a significantly positive correlation between the temporal 
stability and relative biomass (Fig. 7b; Table 2). 
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Fig. 7 (a) Temporal stability of aboveground biomass on non-leguminous species (NL), leguminous species 
(LE), and main leguminous plant species (CS, Caragana sinica; AG, Astragalus galactites; MR, Melissilus 
ruthenicus; OM, Oxytropis myriophylla; TI, Thermopsis lanceolate; AA, Astragalus adsurgens). (b) Linear 
regression between temporal stability and relative biomass. 


Table 2 Correlation coefficients of relative biomass among six leguminous species under mowing (light gray 
shading) and control (dark grey shading) 


Plant species Cs AG MR OM TI AA 
CS - —0.22 0.11 -0.41" -0.69™" —0.21 
AG 0.26 0.58" —0.28 zy 


MR 


Note: CS, Caragana sinica; AG, Astragalus galactites; MR, Melissilus ruthenicus; OM, Oxytropis myriophylla; TI, Thermopsis 
lanceolate; AA, Astragalus adsurgens; *™*, P<0.01 level; *, P<0.05 level. - means no value. 


3.4 TN 


During the initial stage (1980s) of this study, there was no significant difference in TN content 
between mowing and control. However, TN content under mowing significantly decreased compared 
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with control in the 2000s (Fig. 8a). Accordingly, the response ratio (the Napierian logarithm of the 
ratio of N under mowing to N under control) in the 2000s was higher than that in the 1980s (Fig. 8b). 
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Fig.8 Mowing effect on soil total nitrogen content (a) and response ratios in the 1980s and 2000s (b). *, P<0.05 
level; ns, P>0.05 level. Bars are standard errors. 


4 Discussion and conclusions 


4.1 Sensitivity of structure and composition of plant functional group to climate change 


Previous studies have demonstrated that climate change (e.g., temperature and precipitation) 
alters the composition of plant functional group in some ecosystems located in cold regions, such 
as Arctic tundra (Walker et al., 2006), alpine meadows and wetlands (Liu et al., 2017). In the 
temperate grasslands, however, plant functional group is rather resistant to increased temperature 
alone (Hoover et al., 2014). Most of these experiments were conducted within 5 a (Grime et al., 
2008). Shi et al. (2015) found that response of community composition to experimental warming 
in the first 7 a of their study was also resistant in the temperate grasslands, while altered 
environments under long-term climate change could induce changes in the community 
composition (Smith et al., 2009). Our results showed that increased temperature resulted in rapid 
decrease in the relative diversity of leguminous species from 1982 to 2011, with a corresponding 
increase in the relative diversity of non-leguminous species. This is consistent with the findings 
obtained by Butterfield (2015) and Butterfield and Munson (2016), who found that ecological 
communities tracked variation in temperature much more closely, intimating strong responses to 
temperature increases. Three possible reasons may partly be responsible for the rapid loss in the 
relative diversity of leguminous species under long-term climate change. First, in the arid or 
semi-arid grasslands, leguminous species have relatively small proportional abundance in the 
community, making a lower competitive ability for water and nutrient resources compared with 
non-leguminous species. Second, warming-induced reduction of interspecific interactions 
between leguminous species and non-leguminous species may reduce the growth responses of 
leguminous species to elevated temperature (Niu and Wan, 2008). Finally, in the water-limited 
grasslands, leguminous species are more sensitive to changes in water availability (Blessing et al., 
2018). Therefore, reduced water availability in the surface soil under long-term warming (Wan et 
al., 2005) may result in less leguminous species. 

Furthermore, 63.4% of variation in the relative biomass decline of leguminous species could be 
indirectly explained by climatic factors (e.g., temperature). Miller et al. (2007) have predicted that 
the advantage of leguminous over non-leguminous species will be lost as a result of accelerated 
soil nitrogen mineralization and availability under climate warming. Our findings of a tendency 
for leguminous relative biomass to reduce and non-leguminous relative biomass to increase under 
warming partially support the prediction. This may be attributable to the spatial variations in plant 
functional traits, such as plant growth rates and variable drivers of plant nitrogen use efficiency 
between leguminous and non-leguminous species. 

In addition, when leguminous species persist in communities, diversity effects are manifest 
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through their influence on vegetation dynamics and abiotic factors (Finn et al., 2013; Suter et al., 
2015). In this study, the diversity effects of leguminous species significantly decreased from 1982 
to 2011, which was partly explained by the negative effect of elevated temperature. The temporal 
stability of leguminous species was significantly lower than that of non-leguminous species. This 
finding is consistent with previous studies that greater plant diversity in the grasslands may 
improve AGB stability in the face of human or abiotic disturbances (Isbell et al., 2009; Polley et 
al., 2013). This reveals that large species pools likely possess one or more stress-tolerant species 
to resist disturbance (e.g., mowing; Tracy and Sanderson, 2004). Stress-tolerant species can 
compensate for the loss of other species during a disturbance, and thus help stabilize biomass 
across time (Bulleri et al., 2016). For example, increasing diversity can also stabilize communities 
thus increasing component species stability (Downing et al., 2014). In our experiment, the 
sensitivity of leguminous species appeared to weaken species interactions between leguminous 
and non-leguminous species, which can diminish the temporal stability of leguminous species. 
This is likely caused by the weakened species interactions failing to promote niche partitioning or 
facilitation (Isbell et al., 2009). 


4.2 Sensitivity of structure and composition of plant functional group to mowing 


Long-term mowing has significant effects on community structure and composition in the 
grasslands (Shi et al., 2015). In our study, an obvious decrease in AGB and species diversity for 
leguminous and non-leguminous species was observed under mowing from 1982 to 2011. However, 
the magnitude of mowing effect on biomass and diversity was modified by plant functional group, 
with a higher change rate in leguminous species than in non-leguminous species. This may result 
from the same underlying mechanism that mowing leads to significant differences in energy 
partitioning between leguminous and non-leguminous species (Shao et al., 2012). This is 
consistent with the result obtained by Liu et al. (2018), who performed a survey on the changes in 
plant diversity, composition and productivity in the semi-arid grasslands, finding that mowing 
significantly reduced above-ground net primary productivity by 21.7%. However, they also found 
that community structure was not altered by a short-term mowing, which was not in line with our 
results. This indicates that plant community composition changes little under a short-term 
mowing, but the long-term mowing has great effects on the composition of vegetation community, 
such as the rapid loss of leguminous species. This discrepancy may be associated with the 
difference in experimental periods. The positive feedback between aboveground vegetation 
communities and soil micro-organisms might be weakened or offset over the course of our 
long-term mowing treatment. This strongly suggests the decline in the potential for 
overcompensation of different plant functional groups in response to continuous mowing (Liu et 
al., 2017). 

Our results suggest that soil TN content decreases under long-term mowing. However, this 
finding does not agree with a previous study conducted by Robson et al. (2007), who proposed 
that positive feedback between grassland management regimes and dominant plant functional 
strategies maintained high nitrogen availability in the subalpine grasslands. In their study, no 
mowing could actually reduce nitrogen availability and nitrogen transformation process, which 
was correlated with a shift in the vegetation community towards more conservative functional 
strategies and greater dominance by grasses (Robson et al., 2007). This difference with our 
findings may be explained by the following phenomena: (1) vegetation community composition 
in the subalpine grassland is shaped by a colder climate, which is different from the climate of the 
temperate grasslands, therefore, this feedback mechanism may be weakened or even disappear; 
and (2) a short-term experiment may be unable to detect the effect of long-term climate variation 
on biogeochemical cycles. More importantly, the long-term mowing may decrease organic matter 
decomposition and plant nutrient availability, causing a decline in total soil nitrogen in semi-arid 
grasslands (Hobbie, 2015; Rumpel et al., 2015). 


4.3 Implications of mowing- and climate-induced vegetation community changes for the 
semi-arid grasslands 


Overall, our results concerning the responses of leguminous and non-leguminous species to the 
long-term mowing and climatic factors have several implications for sustainable management and 
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utilization strategies in the semi-arid grasslands. First, mowing in the entire grasslands may 
decrease spatial heterogeneity and therefore not be the most appropriate option for maximizing 
plant biomass and species diversity in the grasslands (Baoyin et al., 2014). Kohler et al. (2005) 
suggested that applying regular mowing in the central portion of a meadow but with unmown or 
less frequently mown strips can be a better management alternative. In the long run, this may also 
be a better mowing strategy for the typical semi-arid grasslands, where a variable mowing regime 
is favorable for vegetation community succession and preserving soil nutrient levels (especially 
by leaving unmown strips). Second, our results show that soil TN decreased along with 
leguminous AGB and species diversity. Leguminous species are characterized by their ability to 
biologically fix N2, which has the benefits of being ecologically benign and reducing the use of 
fossil fuels in agriculture. Accordingly, planting leguminous species is usually a more effective 
agronomic practice for supplying nitrogen for agricultural production relative to artificial nitrogen 
application (Zahran, 1999; Jensen and Hauggaard-Nielsen, 2003). We therefore suggest that 
reseeding appropriate leguminous plants into the semi-arid grasslands when mowing may be a 
better strategy to improve nitrogen levels of grassland ecosystems and to thus maintain ecosystem 
biodiversity. 
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Fig.S1 Characteristics of precipitation (a) and temperature (b) during the growing season from 1982 to 2011 
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Fig. S? AGB (a and b) and H (c and d) in non-leguminous species from 1982 to 2011. Inset bar graphs in the 
four panels compare the variable changes between control (black) and mowing (white) from 1982 to 2011. **, 
P<0.01 level; *, P<0.05 level; ns, P>0.05 level. Bars are standard errors. 
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Fig. S3 Biomass/diversity ratios of both leguminous (a) and non-leguminous (b) species under mowing and 
control. **, P<0.01 level. Bars are standard errors. 


